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Abstract: We introduce a new method of delta hedging. In many cases, this 
method results in a lower cost than the Black-Scholes method. To calculate the 
cost of hedging, we develop a Mathematica program. 
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1. Introduction 

We introduce a pricing method of distribution (see Hirashita [4] ) under the condition 
that the risk-free interest rate is equal to the growth rate. This pricing method 
yields a new method of delta hedging. In this paper, we re-calculate several 
examples given in Hull [5] and Luenberger [7], and show that in many cases, the new 
method of hedging results a lower cost than the Black-Scholes method. In order 
to calculate the cost of hedging under this new method, we utilize a Mathematica 
program that include the two-dimensional Newton-Raphson method. It should be 
noted that Mathematica is a programming language. 

2. Well-known results: Delta and gamma with the Black-Scholes 
formula 

We assume that the stock price Y = Se x+rT is lognormally distributed with 
volatility aVT ', where S is the stock price factor, r is the continuously compounded 
interest rate, K is the exercise price of the call option, and T is the exercise period. 
Then, the European call option is given by 
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dF(x) := p(x)dx := e 2^t dx, 
and x £ I :— (— oo, +00). The expectation E of this option is calculated as 
E = / a(x)dF(x) = / (Se x+rT - K)e dx 
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where N(x) = e~ x I 2 /V2n dx is the cumulative standard normal distribution 
function. We set 
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then, the equation E/u — e rT yields the price 

u = SN{d 1 )- Ke- rT N{d 2 ), (2.1) 
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which is the Black-Scholes formula for a European call option (see Hull [5]). Delta 
and gamma are given as follows: 
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delta: A := ^= =JV(di). (2.2) 
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gamma : £ := = e^/C^V^T). (2.3) 



3. New results: Delta and gamma with the simultaneous equations 

It should be noted that inf x£ j a(x) = and J f a{x)dF(x) < oo. As J a r x ^ =0 dF(x) > 
0, the price u and the optimal proportion of investment t u are determined by the 
simultaneous equations 
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(see Corollary 5.1 and Section 6 in Hirashita [4]). Set /3(x) := a(x)t u — ut u + u; 
then, we obtain 

f /j log 0(x)dF(x) =rT + log u, 

As a(x), t„, and w are functions with respect to 5, from 9 (J 7 log/3(x)(ii ;l (x) — logu) /<95 
= 0, we have 
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this implies that 



where 
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Using the well-known formula 

— / F(S, x)dx = F(S, $(S))1/(S)-F(S, <p(S))<p'(S)+ / — (5, x)dx, 



we have 
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e x p(x)dx. 
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From d (j) l/0(x)dF(x) - l/u) /dS = 0, we obtain 

dt u 
dS 

(1 - t u )uWe^ f, jfadFjx) + t u e^ £ gf _ rT {(Se , +rT _^ uU+uf dF(x) - ^ 

-ir-^ydFix) 

Therefore, we can calculate 

gamma: V := = ^(uWe rT ) = uW 2 e 2rT + ue rT ^. (3.3) 

4. Comparison between two delta hedging methods 

When S = 49, K = 50, a = 0.2, r = 0.05, and T = 20/52 (see Section 14.1 in Hull 
[5]), we have 

a(x) = max(49e 1/5 V - 50, 0), dF{x) = J — e ^^dx, 

and the price u = 1.774 with respect to Equation 3.1. At this price, if investors 
continue to invest t u = 0.115 of their current capital, they can maximize the limit 
expectation of the growth rate to e rT = e 1 / 52 = 1.019. Moreover, we have A = 0.448 
and T = 0.0668. 

Meanwhile, the Black-Scholes formula yields u = 2.401, A = 0.522, and T = 
0.0655. 

In the following examples, we consider the cost of hedging for 100, 000 stocks. 

Example 4.1. Table 14.2 in Hull [5] presents a simulation of the delta hedging 
of a sequence of weekly stock prices Si := {49.00, 48.12, 47.37, 50.25, 51.75, 53.12, 
53.00, 51.87, 51.38, 53.00, 49.88, 48.50, 49.88, 50.37, 52.13, 51.88, 52.87, 54.87, 54.62, 
55.87, 57.25} with K = 50, a = 0.2, and r = 0.05. For Si, the cost of hedging is 
$287, 500 according to Equation 3.2, which is 9.2% higher than the cost of hedging 
using the Black-Scholes method ($263,300). 

On the other hand, when K = 65, the cost of hedging is $2, 600 according to 
Equation 3.2, which is 46.9% lower than the cost of hedging using the Black-Scholes 
method ($4,900). 

In addition, when K = 35, the difference between these costs is within 0.1%. 

Example 4.2. Table 14.3 in Hull [5] shows a simulation of the delta hedging 
of a sequence of weekly stock prices S 2 := {49.00, 49.75, 52.00, 50.00, 48.38, 48.25, 
48.75, 49.63, 48.25, 48.25, 51.12, 51.50, 49.88, 49.88, 48.75, 47.50, 48.00, 46.25, 
48.13, 46.63, 48.12} with K = 50, a = 0.2, and r = 0.05. For S 2 , the cost of 
hedging is $247, 900 according to Equation 3.2, which is 3.4% lower than the cost 
of hedging using the Black-Scholes method ($256, 600). 

On the other hand, when K = 65, the cost of hedging is $3, 100 according to 
Equation 3.2, which is 48.3% lower than the cost of hedging using the Black-Scholes 
method ($6,000). 

In addition, when K = 35, the difference between these costs is within 0.1%. 



Example 4.3. We select Table 13.1 in Luenberger [7] of weekly stock prices 
such that S 3 := {35.50, 34.63, 33.75, 34.75, 33.75, 33.00, 33.88, 34.50, 33.75, 34.75, 
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34.38, 35.13, 36.00, 37.00, 36.88, 38.75, 37.88, 38.00, 38.63, 38.50, 37.50} with K = 
35, a = 0.18, and r = 0.1. For S3, the cost of hedging is $274,900 according to 
Equation 3.2, which is 1.3% higher than the cost of hedging using the Black-Scholcs 
method ($271,300). 

On the other hand, when K = 45, the cost of hedging is $1, 300 according to 
Equation 3.2, which is 82.7% lower than the cost of hedging using the Black-Scholes 
method ($7,500). 

In addition, when K = 25, the difference between these costs is within 0.1%. 
5. Program 

In order to confirm that the cost of hedging is indeed $287, 500 in Example 4.1, it is 
necessary to run the following Mathematica program that includes the Newton- Raphson 
method. In order to confirm the other costs of hedging with respect to Equation 
3.2, substitute the Stock, K, sigma, and r data with those given below. 

Stock = {49, 48.12, 47.37, 50.25, 51.75, 53.12, 53, 51.87, 51.38, 
53, 49.88, 48.50, 49.88, 50.37, 52.13, 51.88, 52.87, 54.87, 54.62, 
55.87, 57.25}; 

K = 50; sigma = 0.2; r = 0.05; 
m = 20; shares = 100000; 

a[x_] := If[x < Log[K/S] - r*T, 0, S*Exp [r*T] *Exp [x] - K] ; 

p[x_] := Exp[-(x + sigma~2*T/2)~2/(2*sigma~2*T)]/Sqrt[2Pi*T*sigma~2] ; 

wr = 50; mr = 10; sd = 1000; M = 10~12; 

Price := Module [{u, t}, 

u = NIntegrate [a[x]p[x] , {x, -Infinity, Infinity}, WorkingPrecision 
-> wr, MaxRecursion -> mr, SingularityDepth -> sd]/2; t=0.5; 
Do [{If [u < 1/M, Break []] ; 

f = Exp[NIntegrate[Log[(a[x]t/u - t + l)]*p[x], {x, -Infinity, 
Infinity}, WorkingPrecision -> wr, MaxRecursion -> mr, 
SingularityDepth -> sd] ] - Exp[r*T]; 

g = NIntegrate [(a [x] - u)/(a[x]t - u*t + u) *p [x] , {x, -Infinity, 
Infinity}, WorkingPrecision -> wr, MaxRecursion -> mr, 
SingularityDepth -> sd] ; 

fu = (f + Exp[r*T] )*NIntegrate[-a[x]t/(a[x]t - u*t + u)/u*p[x], 
{x, -Infinity, Infinity}, WorkingPrecision -> wr, MaxRecursion -> 
mr, SingularityDepth -> sd] ; 
ft = 0; 

gu = -NIntegrate [a [x] / (a [x] t - u*t + u)~2*p[x], {x, -Infinity, 
Infinity}, WorkingPrecision -> wr, MaxRecursion -> mr, 
SingularityDepth -> sd] ; 

gt = -NIntegrate [ (a [x] - u)~2/(a[x]t - u*t + u)~2*p[x], {x, 
-Infinity, Infinity}, WorkingPrecision -> wr, MaxRecursion -> mr, 
SingularityDepth -> sd] ; 

ans = Solve [Re [fu] *a + Re[ft]*b == -Re[f] kk Re [gu] *a + Re [gt] *b 
== -Re[g], {a, b}] ; If [ans=={}, Break []] ;a0 = a /. ans[[l]]; 
bO = b /. ans[[l]]; If[Abs[a0] + Abs [bO] < 1/M, Break[]]; u = u + aO; 
t = t + bO; If [u < 0, u = (u - a0)/2]; If [t < 0, t = (t - b0)/2]; 
If [t >= 1, t = ((t - b0)+l)/2] ; 
}, {j2, 1, 100}] ; tO = t; 



5 



u]; 

Print ["{Week, Stock price, Delta, Shares, Cost of shares, 
Cumulative cost, Interest}"]; 
s = 0; cost = 0; interest = 0; 

Do[S = Stock [[j + 1]]; T = (m - j)/52; If [T > 0, u = Price; t = tO, 
u = If [S > K, S - K, 0] ; t =.] ; 

If [T > 0, W = NIntegrate[Exp[x]/((S*Exp[r* T] *Exp [x] - K)t - 
u*t + u)*p[x], {x, Log[K/S] - r*T, Infinity}, WorkingPrecision -> 
wr, MaxRecursion -> mr, SingularityDepth -> sd] , If [ S > K, 
W = 1/(S - K) , W = 0]] ; 

delta = u*Exp[r*T] *W; delta = Round [1000*delta] /1000 . ; 
s2 = shares*delta - s; s = shares*delta; 
cost = cost + Round [s2*S/100] /10 + interest; 
interest = Round [cost*r*l/52*10] /10 . ; 

Print[{j, S, u, delta, s2, Round [s2*S/100] /10 . , Round [cost* 1000] , 
interest}] ; , 
■Cj , 0, m}] ; 

cost = cost - delta*shares*Min[K, Stock [ [m + 1]]]/1000 + (1 - delta)* 
shares*Max [Stock [[m +1]] - K, 0]/1000; 

Print["K = ", K, ", Cost of Hedging =", Round [cost* 1000] ] ; 
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